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Abstract Cosmic ray (CR) transport and acceleration is determined by the proper- 
ties of magnetic turbulence. Recent advances in MHD turbulence call for revisions 
in the paradigm of cosmic ray transport. We use the models of magnetohydrody- 
namic turbulence that were tested in numerical simulation, in which turbulence is 
injected at large scale and cascades to to small scales. We shall address the issue of 
the transport of CRs, both parallel and perpendicular to the magnetic field. We shall 
demonstrate compressible fast modes are dominant cosmic ray scatterer from both 
quasilinear and nonlinear theories. We shall also show that the self-generated wave 
growth by CRs are constrained by preexisting turbulence and discuss the process in 
detail in the context of shock acceleration at supernova remnants and their implica- 
tions. In addition, we shall dwell on the nonlinear growth of kinetic gyroresonance 
instability of cosmic rays induced by large scale compressible turbulence. This gy- 
roresonance of cosmic rays on turbulence is demonstrated an important scattering 
mechanism in addition to direct interaction with the compressible turbulence. The 
feedback of the instability on large scale turbulence should be included in future 
simulations. 



1 Introduction 

The propagation and acceleration of cosmic rays (CRs) is governed by their inter- 
actions with magnetic fields. Astrophysical magnetic fields are turbulent and, there- 
fore, the resonant and non-resonant (e.g. transient time damping, or TTD) inter- 
action of cosmic rays with MHD turbule nce is the accepted principal mechanism 



to scatter and isotropize cosmic rays (see ISchlickeiseri 120021) . In addition, efficient 



scattering is essential for the acceleration of cosmic rays. For instance, scattering of 
cosmic rays back into the shock is a vital component of the first order Fermi accel- 
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eration (see lLongain 1 1 9971) . At the same time, stochastic acceleration by turbulence 
is entirely based on scattering. The dynamics of cosmic rays in MHD turbulence 
holds the key to all high energy astrophysics and related problems. 

We live in an exciting era when we are starting to test fundamental processes 
taking place at the beginning of the Universe, at the event horizon of black holes, 
when the nature of dark matter and dark energy is being probed etc. In the mean 
time, with the launching of the new facilities like Fermi, we have much more obser- 
vational data available than ever before. Using computers many researchers make 
sophisticated complex models to confront the observations in unprecedented details. 
This makes it urgent that we understand the key physical processes underlying as- 
trophysical phenomena, can parameterize them and, if necessary, use as a subgrid 
input in our computer models. 

At present, the propagation of the CRs is an advanced theory, which makes use 
both of analytical studies and numerical simulations. However, these advances have 
been done within the turbulence paradigm which is being changed by the current 
research in the field. Instead of the empirical 2D+slab model of turbulence, numer- 
ical simulations suggest anisotropic Alfvenic modes (an analog of 2D, but not an 
exact one, a s the anisotropy change s with the scale involved) + fast modes or/and 
slab modes ( Cho & LazarianL[2002h . 

At the same time, one should not disregard the possibilities of generation of 
additional perturbations on small scales by CR themselves. For instance , the slab 
Alfvenic perturbation can be created, e.g., via streaming instability (see WentzelL 



1974 ICesarskvl Il980h . These perturbations are present for a range of CRs ener- 



gies (e.g., ~< lOOGeV in interstellar medium for the str eaming instability) ow- 
ing to non-linear damping arisi ng from ambient turbulence (I Yan & Lazarianl 12002 , 



2004 henceforth YL02, YL04, Farmer & Goldreich[ |20Q4|). Instabilities induced by 



anisotropic dis t ribution of CRs were also suggested as a possibility to scatter CRs 
( Lerchel 1967 : Melrose , 1974 ). Gyroresonance instability induced by large scale 
turbulence compressions, can be a n important feedback processe s that creates slab 
mode s to efficiently scatter CRs dLazarian & Beresnvaki l2006t I Yan & Lazarian , 



2011). 



Propagation of CRs perpendicular to mean magnetic field is another important 
problem for which one needs to take into account both large and small scale in- 
teractions in tested models of turbulence. Indeed, if one takes onl y the di f fusion 
along the magnetic fie l d line and fie l d line random walk (FLRW iJokipii , 1 19661 



Jokipii & Parked Il969t iForman et al.L 11974b . compound (or subdiffusion) would 



arise. Whether t he subdiffusion i s real i stic in fact depends on the models of tur- 
bulence chosen (I Yan & Lazarianl 120081 120121) . In this paper we review current un- 
derstandings to this question within the domain of numerically tested models of 
MHD turbulence. 

In what follows, we introduce the basic mechanisms for the interactions between 
particles and turbulence in §2. We discuss the cosmic ray transport in large scale 
turbulence, including both analytical and numerical studies in §3. Applications to 
cosmic ray propagation is presented in §4. In §5, we consider the perpendicular 
transport of cosmic rays on both large and small scales. We shall also discuss the 
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issue of super-diffusion and the applicability of sub-diffusion. In §6, we concentrate 
on the issue of self-confinement in the presence of preexisting turbulence and dwell 
on, in particular, the streaming instability at supernova remnant shocks and its impli- 
cation for CR acceleration. §7, we address the issue of gyroresonance instability of 
CRs and its feedback on large scale compressible turbulence. Summary is provided 
in §8. 



2 Interactions between turbulence and particles 



Basically there are two types of resonant interactions: gyroresonance acceleration 
and transit acceleration (henceforth TTD). The resonant condition is (O—knVfJ. = nQ 
(n = 0,±1,2...), where (0 is the wave frequency, Q = Qo/yis the gyrofrequency of 
relativistic particle, jj. = cos % , where £ is the pitch angle of particles. TTD formally 
corresponds to « = and it requires compressible perturbations. 

The Fokker-Planck equation is generally used to describe the evolvement of the 
gyrophase-averaged distribution function /, 
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where p is the particle momentum. The Fokker-Planck coefficients Du,n,Dnp,D p p 
are the fundamental physical parameters for measuring the s tochastic interactions, 
which are determined by the electromagnetic fluctuations (see lSchlickeiser & Millet , 
19981) : 

Gyroresonance happens when the Doppler shifted wave frequency matches the 
Larmor frequency of a parti cle. In quasi-linear theory (Q LT), the Fokker-Planck 
coefficients are given by (see lSchlickeiser & Milled. 1 19981) 
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where £ = 1 for Alfven modes and £ = k\\/k for fast modes, k m j„ = L~ l , k c — Ho /v f /, 
corresponds to the dissipation scale, m = ymij is the relativistic mass of the pro- 
ton, v± is the particle's velocity component perpendicular to Bo, <p = arctan(k y / k x ) , 
Ji?,& = (jc± ry) / a/2 represent left and right hand polarization. Mij and Ky are the 
correlation tensors of magnetic and velocity fluctuations. 
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From the resonance condition, we know that the most important interaction oc- 
curs at = k^„, s = £2/v\\. This is generally true except for small fi (or scattering 
near 90°). 

TTD happens due to the resonant interaction with parallel magnetic mirror force. 
Particles can be accelerated by when they are in phase with the waves either by 
interacting with oscillating parallel electric field (Landau damping), or by mov- 
ing magnetic mirrors (TTD). When particles are trapped by moving in the same 
speed with waves, an appreciable amount of interactions can occur between waves 
and particles. Since head-on collisions are more frequent than that trailing colli- 
sions, particles gain energies. Different from gyroresonance, the resonance func- 
tion of TTD is broadened even for CRs with small pitch angles. The formal reso- 
nance peak k\\/k = Vpj,/v\\ favors quasi-perpendicular modes. However, these quasi- 
perpendicular modes cannot form an effective mirror to confine CRs because the 
gradient of magnetic perturbations along the mean field direction V||B is small. 
As we will show later in £]3.3I the resonance is broadened in nonlinear theory (see 



Yan & Lazarian, 2008). 



3 Scattering of cosmic rays 

3.1 Scattering by Alfvenic turbulence 



As we discussed in §2, Alfven modes are anisotropic, eddies are elongated along the 
magnetic field, i.e., k± > k^. The scattering of CRs by Alfven modes is suppressed 
first because most turbulent energy goes to k± due to the anisotropy of the Alfvenic 
turbulence so that there is much less energy left in the resonance point k\\ res = 
£2/v\\ ~ r£ l . Furthermore, k± S> k\\ means k± ^> rj) so that cosmic ray particles 
have to be interacting with lots of eddies in one gyro period. This random walk 
substantially decreases the scattering efficiency. The scattering by Alfven modes 
was studied in YL02. In case that the pitch angle % not close to 0, the analytical 
result is 
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where F[a,z] is the incomplete gamma function. The presence of this gamma func- 
tion in our solut ion makes our results orders of magnitude larger than thoseQ in 
Chandran d2000h . who employed GS95 ideas of anisotropy, but lacked the quantita- 
tive description of the eddies. However, the scattering frequency, 



(3) 



1 The comparison was done with the resonant term in Chandran (2000) as the nonresonant term is 
spurious 
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are nearly 10 times lower than the estimates for isotropic and slab model (see 
Fig. HJ/e/f)- It is clear that for most interstellar circumstances, the scattering by 
Alfvenic turbulence is suppressed. As the anisotropy of the Alfven modes is in- 
creasing with the decrease of scales, the interaction with Alfven modes becomes 
more efficient for higher energy cosmic rays. When the Larmor radius of the parti- 
cle becomes comparable to the injection scale, which is likely to be true in the shock 
region as well as for very high energy cosmic rays in diffuse ISM, Alfven modes 
get important. 



3.2 Cosmic ray scattering by compressible MHD turbulence 



As we mentioned earlier, numerical simulations of MHD turbulence supported the 
GS95 model of turbulence, which does not have the "slab" Alfvenic modes that 
produced most of the scattering in the earlier models of CR propagation. Can the 
turbulence that does not appeal to CRs back-reaction (see §4) produce efficient scat- 
tering? 



In the models of ISM turbulence dArmstrong et al.l 1995; Mc Kee & Ostrikei 



2007), where the injection happens at large scale, fast m odes were identified as 



a sca ttering agent for cosmic rays in interstellar medium (lYan & Lazarianl 12002 



2004. These work s made use of the quantitative description of turbulence obtained 
in lCho & Lazarianl d2002h to calculate the scattering rate of cosmi c rays. 



Different from Alfven and slow modes, fast modes are isotropic ( Cho & Lazarianl 



20021) . In deed they are subject to both collisional and collisionless damping. The 
studies in lYan & Lazarianl (120021 12004ft demonstrated, nevertheless, that the scat- 



tering by fast modes dominates in most cases in spite of the damping (see FigQ] 
right). More recent studies of cosmic ray propagati on and acceleration t hat ex - 
plicitly appeal to the effec t of the fast modes inclu de Cassano & Brunettil ( 2005b : 
Brunetti & Lazarianl d2007b ; lYan & Lazarianl d2008b ; lYan et al.l d2008l) . Incidentally, 
fast modes have bee n also identified as pr i mary agents for the acceleration of 
charged dust particles lYan & Lazarianl (120031) : lYan et al.l 020041) . 



3.3 Nonlinear theory of diffusion 

While QLT allows easily to treat the CR dynamics in a local magnetic field system 
of reference, a key assumption in QLT, that the particle's orbit is unperturbed, makes 
one wonder about the limitations of the approximation. Indeed, while QLT provides 



2 On the basis of weak turbulence theory, |Chandran ( 20Q5j) has argued that high-frequency fast 
waves, which move mostly parallel to magnetic field, generate Alfven waves also moving mostly 
parallel to magnetic field. We expect that the scattering by thus generated Alfven modes to be 
similar to the scatter ing by the fast modes cr eated by them. Therefore we expect that the simplified 
approach adopted in lYan & Xa zarian (2004) and the papers that followed to hold. 
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Fig. 1 Left: rate of CR scattering by Alfven waves versus CR energy. The lines at the top of 
the figure are the acc epted estimates obtained for Kolmogorov turbulence. The dotted curve is 
from Chandran (2000). The analytical calculations are given by the solid line with our numerical 
calculations given by crosses; Right: the scattering by fast modes, dashed line represents the case 
without damping for fast modes included, the solid and dash-dot line are the results taking into 
account collisionless damping. 



simple physical insights into scatterin g, it is known to have problems. For in s tance, 
it fails in trea t ing 90° scattering (see Volkll 19731: IVoelkl 1 19751: Ijones et all Il973 , 
19781: lOwensi 1 1974 iGoldsteini Il97fi iFelice & Kulsru dl200lb and perpendicular 



transport (see lKota & Jokipiil bOOOl: iMatthaeus et all 12003b 

Indeed, many attempts have been ma de to improve the QLT and various non- 
linear theories have been attempted (see lDupreelll966l Volk 1973, 1975, Jones, 
Kaiser & Birmingham 1973, Goldstein 1976). Currently we observe a surge of inter- 
est in findi ng way to go beyond QL T. Examples include the nonlinear guiding center 
theory (seejMatthaeus et aUl2003l). second-order quasilinear theory ( Shalchi , 2005 
!ebbl" ' 



Qin, 2007; le Roux & Webb, 2007), etc. Most of the analysis were limited to tradi- 



tional 2D+slab models of MHD turbulence. An important step was taken in Yan & 
Lazarian (2008), where non-linear effect was accounted for in treating CR scatter- 
ing in the type of MHD turbulence that are supported by numerical simulations. The 
results have been a pplied to both solar flares (Yan, Lazarian & Petrosian 2008) and 



grain acceleration (IHoang et al.l 120121) . Below, we introduce the nonlinear theory 



and their applications to both particle transport and acceleration in incompressible 
and compressible turbulence based on the results from Yan & Lazarian (2008). 

The basic assumption of the quasi-linear theory is that particles follow unper- 
turbed orbits. In reality, particle's pitch angle varies gradually with the variation 
of the magnetic field due to conservation o f adiabatic invariant v\/B , where B is 
the total strength of the magnetic field (see Landau & Lifshitzl 1975b . Since B is 
varying in turbulent field, so are the projections of the particle speed Vj_ and V||. 
This results in broadening of the resonance. The variation of the velocity is mainly 
caused by the magnetic perturbation &B| in the parallel direction. This is true even 
for the incompressible turbulence we discussion in this section. For the incompress- 
ible turbulence, the parallel perturbation arises from the pseudo-Alfven modes. The 
perpendicular perturbation SB± is higher order effect, which we shall neglect here. 
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The propagation of a CR can be described as a combination of a motion of its 
guiding center and CR's motion about its guiding center. Because of the dispersion 
of the pitch angle Afi and therefore of the parallel speed 4v||, the guiding center is 
perturbed about the mean position < z >= vjlt as they move along the field lines. 
As a result, the perturbation 8B(x,t) that the CRs view when moving along the 
field gets a different time dependence. The characteristic phase function e**ll z w of 
the perturbation 8B(x,t) deviates from that for plane waves. Assuming the guiding 
center has a Gaussian distribution along the field line, 

, (:-<z>) 2 

m = -=- e^r, (4) 

V27T cr- 
one gets by integrating over z, 

dze ik W z f(z)=e lk \\ < ' > e- k i a ' /2 . (5) 



The first adiabatic invariant gives us 




(6) 



Insert the Eq.© into the expression of (see Volk 1975. lYan & Lazarian 



20041) . we obtain 



[ d 3 kY^R n {k\\v\\- co±nQ) 
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Following are the definitions of the parameters in the above equation. £2 , ji are the 
Larmor frequency and pitch angle cosine of the CRs. J„ represents Bessel function, 
and w — kj_v±/Q = k±LR^l — jx 2 , where R = v/(Ql) is the dimensionless rigidity 
of the CRs, L is the injection scale of the turbulence. kx,k\\ are the components of the 
wave vector k perpendicular and parallel to the mean magnetic field, co is the wave 
frequency. / (k) is the energy spectrum of the Alfven modes and 7 M (k) represents 
the energy spectrum of magnetosonic modes. In QLT, the resonance function R n = 
n8(k^ V|| — (OzknQ). Now due to the perturbation of the orbit, it should be 

R„(k v -co±nQ) 

= [™ dte (k ^ +na ~ m), ~^ k l <Av r >t2 
Jo 
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where Ma = SV/va = SB/Bq is the Alfvenic Mach number and va is the Alfven 
speed. We stress that Eqs. (I7l8l l are generic, and applicable to both incompressible 
and compressible medium. 

For gyroresonance (n = ±1,2,...), the result is similar to that from QLT for 
H 3> A/1 = 4v||/v. In this limit, Eq.(|7]i represents a sharp resonance and be- 
comes equivalent to a (5-function when put into Eq.(|7]i. In general, the result is 
different from that of QLT, especially ata-> 90°, the resonance peak happens at 
foi res ~ Q/Av in contrast to the QLT result kn res ~ &/v\\ —> °° ■ We shall show be- 
low, that due to the anisotropy, the scattering coefficient is still very small if 
the Alfven and the pseudo-Alfven modes are concerned. 

On the other hand, the dispersion of the vii means that CRs with a much wider 
range of pitch angle can be scattered by the pseudo-Alfven modes through TTD (n = 
0), which is marginally affected by the anisotropy and much more efficient than the 
gyroresonance. In QLT, the projected particle speed is comparable to phase speed 
of the magnetic field compression. This means that only particles with a specific 
pitch angle can be scattered. For the rest of the pitch angles, the interaction is still 
dominated by gyroresonance, which efficiency is negligibly small for the Alfvenic 
anisotropic turbulence (see £13 - lb - With the resonance broadening, however, wider 
range of pitch angle can be scattered through TTD, including 90°. 



3.4 Results from test particle simulations 



We live in an era when we can test various processes in astrophysics and numerical 
studies have become an important part of theoretical efforts. Test particle simula - 
tion has been used to study CR scattering and transport lGiacalone & Jokipiil (119991) : 
Mace et al. (120001) . The aforementioned studies, however, used synthetic data for 
turbulent fields, which have several disadvantages. Creating synthetic turbulence 
data which has scale-dependent aniso tro py with respect to the local magnetic field 
(as observed in lCho & Vishniad2OO0l and Maron & GoldreichkobH) is difficult and 
has not been realized yet. Also, synthetic data normally uses Gaussian statistics and 
delta- correlated fields, which is hardly ap propriate for description of strong turbu- 
lence. Beresnyak. Yan. & Lazarianl ( 201 1 ) for the first time used the results of direct 
numerical MHD simulations a s the input data for parti cle scattering simulations. 

Both the results obtained in lBeresnvak et al.l (1201 II) and recent ones (Xu & Yan, 
in preparation) are quite encouraging with good correspondence with the analytical 
predictions. We briefly summarize the results here. As shown in Figf2] particles 's 
motion is diffusive both along the magnetic field (x direction) and across the field 
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(y direction). Moreover, the scattering coefficient is almost constant with energy as 
revealed in Fig.\2\right. This is exac tly due to the nonlinear TTD interaction, which 
is scale independent, as predicted in lYan & Lazarianl (120081) . 




Fig. 2 Left: Diffusive behavior of the particles in the tracing simulations. For different rt we show 
ensemble-averaged square deviations, which are proportional to time, x and y are measured in units 
of the cube size; right: dimensionless CR scattering coefficient D^^/Q for the case of [l = 0.71 
vs Larmor radius expressed in cube size units (solid line). It is dominated by magnetic bottles 
formed by slow mode, this is why D^^ (dimensional scattering frequency) is almost constant, i.e. 
independent on particle's energy. For comparison, we plot various theoretical predicti ons: QLT 
prediction for Alfven and slow mode (dashed) (from Beresnvak, Yan, & Lazarian, 201 1). 



4 Cosmic ray propagation in Galaxy 

The scattering by fast modes is influenced by the medium properties as the fast 
modes are subject to linear damping, e.g., Landau damping. Using the approach 
above we revisit the problem of the CR propagation in the selected phases of the ISM 
(see Table |4] for a list of fiducial parameters appropriate for the idealized phasefl) 
assuming that the turbulence is injected on large scales. 



4.1 Halo 



In Galactic halo (see Table @), the Coulomb collisional mean free path is ~ lOpc, 
the plasma is thus in a c ollisionless regime. The cascading rate of the fast modes is 
dCho & Lazarianl 12002b. 



3 The parameters of idealized interstellar phases are a subject of deba te. Recently, even the entire 
concept of the phase being stable entities has been challenged (see iGazol et"aU I2007L and ref. 
ther ein). Indeed differen t parts of interstellar medium can exhibit variations of these parameters 
(see lWolfire et al.l l2003. and ref. therein) 
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ISM 


halo 


HIM 


WIM 


WNM 


CNM 


DC 


T(K) 


2x 10 6 


1 x 10 6 


8000 


6000 


100 


15 


c.s(km/s) 


130 


91 


8.1 


7 


0.91 


0.35 


n(cm~ 3 ) 


10- 3 


4x 10~ 3 


0.1 


0.4 


30 


200 


l,„f p (cm) 


4x 10 19 


2x 10 18 


6x 10 12 


8x 10 11 


3x 10 6 


10 4 


L(pc) 


100 


100 


50 


50 


50 


50 


BQiG) 


5 


2 


5 


5 


5 


15 


P 


0.28 


3.5 


0.11 


0.33 


0.42 


0.046 



damping collisionless collisional collisional neutral-ion neutral-ion neutral-ion 



Table 1 The parameters of idealized ISM phases and relevant damping. The dominant damping 
mechanism for turbulence is given in the last line. HIM=hot ionized medium, CNM=cold neutral 
medium, WNM=warm neutral medium, WIM=warm ionized medium, DC=dark cloud. 




2 - WIM collisionless 

1 °0 20 40 60 80 ? < , ' 10 ! 10* 10 5 

E t (GeV) 

Fig. 3 Left: The turbulence truncation scales in Galactic halo and warm ionized medium 
(WIM). The damping curves flattens around 90° due to field line wandering (dotted lines, see 
I Yan & Laza rian 2004; Lazarian. Vishniac, & Chol2004l) : For WIM, both viscous and collisionless 
damping are applicable; right: The mean free paths in two different phases of ISM: halo (solid 
line) and WIM (dashed line). At lower energi es (~< lOOGeV ), the different dependence in WIM 
is owing to the viscous damping (from lYan & L azarian. 2 008I) , 



X^ = (k/L) l / 2 SV 2 /V ph 
By equating it with the collisionless damping rate 
/nfisin 2 9 



2 cos 9 
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we obtain the turbulence truncation scale k c : 

■ exp 



4Mjm i cos 2 

k c L~ 



7tm e p sin 9 



2m? 



j8to, cos 2 9 



J3 cos 2 9 



(9) 



(10) 



(11) 



where j3 = P gas /Pmag- 

The scale k c depends on the wave pitch angle 9, which makes the damp- 
ing anisotropic. As the turbulence undergoes turbulent cascade or/and the waves 
propagate in a turbulent medium, the angle 9 is changing. As discussed in YL04 
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the field wandering defines the spread of angles. During one cascading time, 
the fast modes propagate a distance vx cas and see an angular deviation tan 89 ~ 

tan 2 +tan 2 89±, which is 



tan 89 



M\ cos 9 



M\ sin 2 9 



1/3 



\ 27 (kL) 1 / 2 \ kL 



(12) 



As evident, the damping scale given by Eq.dTTb varies considerably especially when 
0^0 and 9 — > 90°. For the quasi-parallel modes, the randomization (°= (fcL) -1 / 4 ) 
is negligible since the turbulence cascade continues to very small scales. On small 
scales, most energy of the fast modes is contained in th ese quasi-parallel modes 
( I Yan & Lazarianl l2004t iPetrosian. Yan. & Lazarian 20061) . 

For the quasi-perpendicular modes, the damping rate (EqJTQb should be averaged 
over the range 90° — 89 to 90°. Equating Eq.® and Eq.dTob averaged over 89, 
we get the averaged damping wave number (see FigJ^/e/f). The field line wander- 
ing has a marginal effect on the gyroresonance, whose interaction with the quasi- 
perpendicular modes is negligible (YL04). However, TTD scattering rates of mod- 
erate energy CRs (< lOTeV) will be decreased owing to the increase of the damping 
around the 90° (see Figj3]/e/0- For higher energy CRs, the influence of damping is 
marginal and so is that of field line wandering. 




lhe l^Ll result in the range jj, > Aji provides a good approximation tor our 
calculations using the non-linear theory (lYan & Lazarian 120081) . For CRs with 
sufficiently small rigidities, the resonant fast modes (k res ~ l/(RjJ.)) are on small 
scales with a quasi-slab structure (see Figj3]/e/f). For the scattering by these quasi- 
parallel modes, the an alytical result that follows from QLT approximation (see 
Yan & Lazarian , 2004 ) for the gyroresonance is0 



4 It can be shown that the QLT result follows from our more general result Eq.(??) if we put 
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(13) 



where tan0 c = fc_L lC /fc|i rejS . This justifies our use of the analytical approximation 
above. 

Once we kn ow the funct ional form of the D^, we can obtain the corresponding 
mean free path dEarll,ll974l) : 



d/i 



v(l-n 



2\2 



MM 



(14) 



where is the co ntribution from TTD in teraction and can be obtained using the 



nonlinear theory (see lYan & Lazarianll2008l and also $3.3b with the inertial range of 
fast modes determined for the local medium (see, e.g.QT|in the case of collisionless 
damping). 

The mean free path is sensitive to the scattering by gyroresonance at small pitch 
angles, due to the influence of damping on the fast modes on small scales. Fig|4] 
shows the pitch angle diffusion of CRs with different energies due to the TTD and 
gyroresonance. 

The weak dependence of the mean free path (see Figf3] right) of the moderate 
energy (e.g< ITeV) CRs in halo results from the fact that gyroresonance changes 
marginally with the CR energy (see Figj4). This is associated with the damping in 
collisionless medium. We expect that similar flat dependence can happen in any 
collisionless medium. This can be a natural explanation of the puzzling "Palmer 
Concensus", the same trend observed in solar wind. 



4.2 Warm Ionized Medium 

In warm ionized medium, the Coulomb collisional mean free path is l m f p = 6 x 
10 12 cm and the plasma /3 ~ 0.11. Suppose that the turbulence energy is injected 
from large scale, then the compressible turbulence is subjected to the viscous damp- 
ing besides the collisionless damping. By equating the viscous damping rate with 
the cascading rate (Eq|9]l, we obtain the following truncation scale, 

k cL =J (15) 

C l(l-3^ 2 )-M»l 



where x c 



6pSV 2 L 
VoV A 



2 



! , Tjo is the longitudinal viscosity. In the low /3 regime, the mo- 



tions are primarily perpendicular to the magnetic field so that dv x /dx = h jn ~ B /B. 
The longitudinal viscosity enters here as the result of distortion of the Maxiwellian 
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distribution (see lBraginskiilll965b . The transverse energy of the ions increases dur- 
ing compression because of the conservation of adiabatic invariant v 2 ± /B. If the rate 
of compression is faster than that of collisions, the ion distribution in the momentum 
space is bound to be distorted from the Maxiwellian isotropic sphere to an oblate 
spheroid with the long axis perpendicular to the magnetic field. As a result, the 
transverse pressure gets greater than the longitudinal pressure, resulting in a stress 
~ T]odv x /dx. The restoration of the equilibrium increases the entropy and causes 
the dissipation of energy. 

The viscous damping scale is compared to collisionless cutoff scale (EqfTTI) in 
Fig.??/e/f. As shown there, both viscous damping and collisionless damping are im- 
portant in WIM. Viscous damping is dominant for small 9 and collisionless damp- 
ing takes over for large 9 except for 9 = 90°. This is because collisionless damping 
increases with 9 much faster than the viscous damping. For sufficiently small wave 
pitch angles, the viscous damping is too small to prevent the fast modes to cascade 
down to scales smaller than the mean free path l m f p . Because of the similar quasi- 
slab structure on small scales, Eq.(IT3b can be also applied in WIM. The results are 
illustrated in FigH] Compared to the case in halo, we see that the qualitative differ- 
ence stands in the gyroresonance. This is because gyroresonance is sensitive to the 
quasi-slab modes whose damping differs in halo and WIM. 



4.3 Other phases 



In hot ionized medium (HIM), the plasma is also in collisionless regime, but the 
density is higher and the plasma beta is larger than 1. The damping by protons 
thus becomes substantial especially at small pitch angles. The damping truncates 
the turbulence at much larger scales than the gyroscales of the CRs of the energy 
range we consider. No gyroresonance can happen and some other mechanisms are 
necessary to prevent CRs streaming freely along the field. The turbulence injected 
from small scales might play an important role (see §5). 

In partially ionized gas one shoul d take into account an additional damping that 
arises from ion-neutral colli sions (see lKulsrud & Pearcetl 19691: Lithwick & Goldreichi 



2001: ILazarian et al.Ll2004l) . In the latter work a viscosity-damped regime of turbu- 



lence was predicted at scales less the scale k~^ mb at which the ordinary magnetic 
turbulence is dam ped by ionic viscosity. The corresponding numerical work, e.g., 
Cho et al.1 ( 12002b testifies that for the viscosity-damped regime the parallel scale 
stays equal to the scale of the ambipolar damping, i.e., k\\ = k cam b, while k± in- 
creases. In that respect, the scattering by such magnetic fluctuations is analogous 
to the scattering induced by the weak turbulence (see §2.3). The difference stems 
from the spectrum of k± is shallower than the spectrum of the weak turbulence. The 
predicted values of th e spectrum for the viscosity-damped turbulence E(k±) ~ k 



-l 



dLazarian et al.L 120041) are in rough agreement with simulations. More detailed stud- 



ies of scattering in partially ionized gas will be necessary. 
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5 Perpendicular transport 



In this section we deal with the diffusion perpendicular to mean magnetic field. 

Propagation of CRs perpendicular to the mean magnetic field is another im- 
portant problem in which QLT encounters serious difficulties. Compound diffu- 
sion, resulting from the convolution of diffusion along the magnetic field line 
and diffusion of field line perpendicular to mean field dire ction, has been in 



yoked to discuss transport of cosmic rays i n the Milky Way (iGetmantsevL 11963 



Lingenfelter. Ramatv. & Fiskl 11971b lAllanl 119721) . The role of compound diffu 



sion in the acceler atio n of CRs at quasi- perpendicular shocks were investigated by 



Duffy et all d 19951) and lKirk et all dl996l) . 



Indeed, the idea of CR transport in the direction perpendicular to t he mean mag- 
netic field being dominated by the field l ine random walk (FLRW, IJokipiil 1 1966 



Jokipii & Parker! 1969 ; Forman et al 1 19741 can be easily justified only in a restricted 



situation where the turbulence perturbations are small and CRs do not scatter back- 
wards to retrace their trajectories. If the latter is not true, the particle motions are 
subdiffusive, i .e., the squared distance diffused growing as not as t but a s t a , CC < 1, 
e.g.. a = 1/2 dKota & JokipiilbOOOtlivlace et alll200d : IOin et alll2002h . If true, this 
could indicate a substantial shift in the paradigm of CR transport, a shift that surely 
dwarfs a modification of magnetic turbulence model from the 2D+slab to a more 
simulation-motivated model that we deal here. 

It was also proposed that with substantial transverse structu re, i.e., transverse 



2002). 



displacement of field lines, perpendicular diffusion is recovered (lOin et al 
Is it the case of the MHD turbulence models we deal with? 

How realistic is the subdiffusion in the presence of turbulence? The answer for 
this question apparently depends on the models of turbulence chosen. 

Compound diffusion happens when particles are restricted to the magnetic field 
lines and p erpendicular transport i s solely due to the random walk of field line wan- 
dering (see iKota & Jokipii , l2000h . In the three-dimensional turbulence, field lines 
are diverging away due to shearing by the Alfven m odes (see lLazarian & Vishniac , 
19991 Naravan & Medvedev, 2001; Lazarian, 2006 . Since the Larmor radii of CRs 
are much larger than the minimum scale of eddies l± m i n , field lines within the CR 
Larmor orbit are effectively diverging away owing to shear by the Alfvenic turbu- 
lence. The cross-field transport thus results from the deviations of field lines at small 
scales, as well as field line random walk at large scale (> min[L/M^,L]). 

Both observation of Galactic CRs and solar wind indicate that the diffusion of 

CRs perpendicular to magneti c field is comparable to parallel diffusion dGiacalone & Jokipii , 
1999tlMaclennan et alll200ll) . Why is that? 



Most recently the diffusion in magneti c fields was considered for thermal parti- 
cles in Lazarian (2006), for cosmic rays in Yan & Lazarian (120081). In what follows 
we present the results based on the studies in I Yan & Lazarianld2008l) . 

For perpendicular diffusion, the important issue is the reference of frame. We 
emphasize that we consider the diffusion perpendicular to the mean field direction 
in the global reference of frame. 
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High Ma turbulence: High Ma turbulence corresponds to the field that is easily 
bended by hydrodynamic motions at the injection scale as the hydro energy at the 
injection scale is much larger than the magnetic energy, i.e. pvl ^> B 2 . In this case 
magnetic field becomes d yna micall y important on a much smaller scale, i.e. the 
scale Ia = L/M^ IEazarian (see 20061) . If the parallel mean free path of CRs hi <C /a, 
the stiffness of B field is negligible so that the perpendicular diffusion coefficient 
is the same as the parallel one, i.e., D± = D« ~ 1/3Aiiv. If X« 3> Ia> the diffusion 
is controlled by the straightness of the field lines, and D± =Z)| w 1 /3/aV, Ma > 
1 , A|| > Ia- The diffusion is isotropic if scales larger than Ia are concerned. 

Low Ma turbulence: In the magnetically dominated case, i.e. the field that cannot 
be easily bended at the turbulence injection scale, individual magnetic field lines 
are aligned with the mean magnetic field. The diffusion in this case is anisotropic. 
If turbulence is injected at scale L it stays weak for the scales larger than LM^ and 
it is strong at smaller scales. Consider first the case of A|| > L. The time of the 
individual step is L/vu, then Dj_ ~ 1 /3LvM*, Ma < 1. This is similar to the case 
discussed in the FLRW model (Jokipii 1966). However, we obtain the dependence 
of M\ instead of their M\ scaling. In the opposite case of A|| < L, the perpendicular 
diffusion coefficient is D±_ w D^M^, which coincides with the result obtained for the 
diffusion of thermal electrons in magnetized plasma (ILazarianl 120061) . This is due 
to the anisotropy of the Alfvenic turbulence. 



5.2 Superdiffusion on small scales 



The diffusion of CR on the scales <C L is different and it is determined by how 
fast field lines are diverging away from eac h other. The mean deviation of a field 
in a d istance 8x is proportional to [5z] 3 / 2 ( Lazarian & Vishniac , 19991: Lazariani 
2006b. same as Ri chardson diffusion in the case of hydro dynamic turbulence (see 
Evinket aUl2011l) . Following the argument, we showed in I Yan & Lazariani d2008l) 
that the cosmic ray perpendicular transport is superdiffusive. The reason is that there 
is no random walk on small scales up to the injection scale of strong MHD turbu- 
lence (LM\ for Ma < 1 and Ia for M a > 1). This can well ex plain the recently 
observed super-diffusion in solar wind (IPerri & Zimbardo l2009h . 



5.3 Is there subdiffusion? 

The diffusion coefficient DnMJ we obtained in the case of Ma < 1, means that 
the transport perpendicular to the dynamically strong magnetic field is a diffu- 
sion, rather than the subdiffusion as discussed in a number of recent papers. This is 
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also supported by test particle simulations (IBeresnvak. Yan. & Lazarianll201 1 



Fig^ight), Let us clarify this point by obtaining the necessary conditions for the 
subdiffusion to take place. 

The major implicit assumption in subdiffusion (or compound diffusion) is that 
the particles trace back their trajectories in x direction on the scale 8z- When is 
it possible to talk about tracing particle trajectories back? In the case of random 
motions at a single scale only, the distance over which the p article trajectories 
get uncorrected is given by the iRechester & Rosenbluthl (119781) model. Assuming 
that the damping scale of the t urbulence is larger that the CR Larmor radius, the 
IRechester & Rosenbluth dl978l) model , when generaliz ed to anisotropic turbulence 
provides dNaravan & Medvedevl l200ll iLazariaru, 120061) Lrr = ^||., m „ln(/x,mm/'"L«r) 
where l\\ m i„ is the parallel scale of the cut-off of turbulent motions, l±. m j n is the 
corresponding perpendicular scale, r^,,. is the CR Larmor radius. The assumption of 
fLar < l±.min can be valid, for instance, for Alfvenic motions in partially ionized gas. 
However, it is easy to see that, even in this case, the corresponding scale is rather 
small and therefore subdiffusion is not applicable for the transport of particles in 
Alfvenic turbulence over scales ^> l\\ >m i n . 

If riar > l_Lmin> as it i s a usual case for Alfven motions in the phase of ISM with 
the ionization larger than w 93%, w here the Alfvenic motions go to the thermal 



partic le gyroradius (see estimates in iLithwick & GoldreichL 1200 It lLazarian et al 



2004), the subdiffusion of CR is not an applicable concept for Alfvenic turbu- 



lence. This does not preclude subdiffusion from taking place in particular models o f 
magnetic perturbations, e.g. in the slab model considered in lKota & Jokipiil feOOd) . 
but we believe in the o mnipresence of Alfvenic turbulence in interstellar gas (see 



Armstrong et all 119951) . 



6 Streaming Instability in the Presence of Turbulence 



When cosmic rays stream at a velocity much larger than Alfven velocity, they can 
excite by gyroresonance MHD modes which in turn scatter cosmic rays back, thus 
increasing the amplitude of the resonant mode. This runaway process is known as 
streaming instability. It was claimed that the insta bility could provi de confinement 



for cosmic rays with energy less than - 10 2 GeV dCesarskyl fl980l) . However, this 



was calculated in an ideal regime, namely, there was no background MHD turbu- 
lence. In other words, it was thought that the self-excited modes would not be ap- 
preciably damped in fully ionized ga^. 

This is not true for turbulent medium, however. Yan & Lazarianl ( 2002 ) pointed 
out that the streaming instabili ty is partially suppres sed in the presence of back- 



ground turbulence (see more in lLazarian et all [2003 ). More recently, detailed cal 



culations of the streami ng instability in the presence of background Alfvenic turbu- 
lence were presented in Farmer & Goldreichl (120041) . The growth rate of the modes 



5 We neglect the nonlinear Landau damping, which is suppressed in turbulence due to decrease of 
mean free path. 
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Table 2 The notation we used in this section 



A 


Normalized wave amplitude SB j Bo 


a 


hardening of the CR spectrum at the shock front 


Bo 


mean magnetic field at the shock in the later Sedov phase 


"CflV 


inercloud magnetic field strength 


SB 


wave amplitude 


c 


light speed 


d 


distance of the molecular cloud from observer 


D 


diffusion coefficient of CRs 


E 


CR energy 


Esn 


supernova explosion 


f 


distribution function of CRs 


fit 


fraction of energy transferred from parent protons to pions 


k 


wave number 


K(t) 


Normalization factor of CR distribution function 


L 


the injection scale of background turbulence 


m 


proton rest mass 


M, 


cloud mass 


n 


intercloud number density 


Ny 


7 ray flux 


P 


CR's momentum 


Pma.x 


the maximum momentum accelerated at the shock front 


PcR 


CR pressure 


q 


charge of the particle 


r 


distance from SNR centre 


Rc 


the distance of the molecular cloud from the SNR centre 


r s 


Larmor radius of CRs 


R d 


diffusion distance of CRs 


Rsh 


shock radius 


Respi tesp 


the escaping distance/time of CRs 


S 


ID spectrum index of CR distribution 


t 


time since supernova explosion 


tage 


the age of SNR 


hed 


the time at which SNR enters the Sedov phase 


U 


shock speed 


Ui 


initial shock velocity 


V 


particle speed 


Vs 


streaming speed of CRs 


w 


wave energy 


a 


power index of D with respect to particle momentum p 


X 


reduction factor of D with respect to Djsm 


8 


power index of p max with respect to t 




fraction of SN energy converted into CRs 




a numerical factor in Eq.?? 


r cr 


the growth rate of streaming instability 




wave damping rate 


K 


ratio of diffusion length to shock radius 


n 


the Larmor frequency of non-relativistic protons 




cross section for pp collision 


k 


the ratio of CR pressure to fluid ram pressure 
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of wave number k is dLongair , 119971) . 



j- ti\ /-> ^) / 1 , v stream\ 

I^ r (^)=X2 (-1+— - — ), (16) 

n p V A 

where N(>E)is the number density of cosmic rays with energy > E which resonate 
with the wave, n p is the number density of charged particles in the medium. The 
number density of cosmic rays near the sun is Af(>£)~2x 10- 10 (£/GeV)- 16 cm 3 sr 
dWentzelL[l974l) . 

Interaction with fast modes was considered by lYan & Lazarianl faOOA Such an 
interaction happens at the rate ~ (k/L)~ l ' 2 Vph/V 2 . By equating it with the growth 
rate Eq.([T6ll. we can find that the streaming instability is only applicable for particles 
with energy less than 

y max ~ 1.5 x 10- 9 [« / 7 1 (y p/l /y)(LvX2o/y 2 )°' 5 ] 1/1 ' 1 , (17) 

which for HIM, provides ~ 20GeV if taking the i njection speed to be V ~ 25 km/s. 



Similar result was obtained with Alfven modes bv lFarmer & Go ldreichl (120041) . 

Magnetic field itself is likely to be amplified through an inverse cascade of mag- 
netic energy at which perturbations created at a particular k diffuse in k space to 
smaller k thus inducing inverse cascade. As the result, the magnetic perturbations 
at smaller k get larger than the regular field. Consequently, even if the instability is 
suppressed for the growth rate given by Eq. ([TBI it gets efficient due to the increase 
of perturbations of magnetic field stemming from the inverse cascade. The precise 
picture of the process depends on yet not completely clear details of the inverse 
cascade of magnetic field. 

Below, we present the application of the current understanding of the interaction 
between the streaming instability and the background turbulence to the modeling 
of the gamma ra y emission from molecular clouds near SNRs (see more details in 



Yan et al.L 120121) . We shall treat the problem in a self-consistent way by comparing 
the streaming level that is allowed by the preexisting turbulence and the required 
diffusion for the CRs. 



6.1 Application to CR acceleration at the shocks 

Diffusive shock acceleration of energetic CR particles relies on the crucial process 
of amplification of MHD turbulence so that particles can be trapped at the shock 
front long enough to be accelerated to the high energy observed. One of the most 
popular scenarios that has been adopted in the literature is the streaming instability 
generated by the accelerated particles. However, in the highly nonlinear regime the 
fluctuations of magnetic field arising from the streaming instability get large and the 
classical treatment of the streaming instability is not applicable. We circumvent the 
problem by proposing that the field amplification we consider does not arise from 
the streaming instability, but is achieved earlier through other processes, e.g. the 
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interaction of the shock precursor with density perturb ations preexisting in the in 



terstellar medium (IBeresnvak. Jones. & Laza rian, 2009). Due to the resonant nature 



of the streaming instability, the perturbations SB arising from it are more efficient in 
scattering CRs compared to t he large scale fluctuatio ns produced by non-resonant 



mechanisms, e.g. the one in IBeresnvak et al.l (120091) . Therefore in this paper, we 



limit our discussions to the regime of SB ~< Bo, where Bq is the magnetic field that 
has already been amplified in the precursor region^. 

When particles reach the maximum energy at a certain time, they escape and the 
growth of the streaming instability stops. Therefore we can obtain the maximum 
energy by considering the stationary state of the evolution. The steady state energy 
density of the turbulence W (k) at the shock is determined by 

(u±v A )vw(k) = 2(r cr -r d )w(k), (is) 

where U is the shock speed, and the term on the l.h.s. represents the advection of 
turbulence by the shock flow, \>a = Bq/ V \nnm and n are the Alfven speed and 
the ionized gas number density of the precursor region, respectively. The plus sign 
represents the forward propagating Alfven waves and the minus sign refers to the 
backward propagating Alfven waves. The terms on the r.h.s. describes the wave am- 
plification by the streaming instability and damping with T c { as the corresponding 
damping rate of the wave. The distribution of accelerated particles at strong shocks 
is f{p) °< p~ 4 . If taking into account the modification of the shock structure by 
the accelerated particles, the CR spectrum becomes harder. Assume the distribu- 
tion of CRs at the shock is f o(p) x p~ 4+a . The nonlinear growth was studied by 
Ptuskin & Zirakashvilil(l2005b . 



The generalized growth rate of streaming instability is 



12^Vv A Vl+A 2 
c 2 k 



/ dpp 




D 


df 


' Pres 






dx 



(19) 



where q is the charge of the particle, c is the light speed, p res = ZeB^X +A 2 /c/k res 
is the momentum of particles that resonate with the waves. A = SB/Bq is wave 
amplitude normalized by the mean magnetic field strength Bq. 

D = Vl+A 2 vr g /3/A 2 (> k res ) (20) 

is the diffusion coefficient of CRs, v and r ? are the velocity and Larmor radius of 
the CRs. In the planar shock approximation, one gets the following growth rate of 
the upstream forward moving wave at x=0, 



6 The eff ective Bp is therefore renor malized and can be much larger than the typical field in ISM 
(see, e.g jDiamond & Malko\i2007l) . 
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r cr {k) 



C cr ^U 2 (U + v A )k 



l-a 



{\+A 2 )^-«)l 2 cv A ${p max )r% 



(21) 



where C cr = 9/2/(4 — a)/ (2 — a), rg = mc 2 /q/Bo, where t, measures the ratio of 
CR pressure at the shock and the upstream momentum flux entering the shock front, 
m is the proton rest mass, and p max is the maximum momentum accelerated at the 
shock front. H(p) is the Heaviside step function. 

The linear damping is negligible since the medium should be highly ionized. 
In fully ionized gas, there is nonlinear Landau damping, which, however, is sup- 
pres sed due to the reduction of particles' mean free path in the turbulent medium 
(see lYan & Lazarianl 1201 ll) . We therefore neglect this pr ocess here. Backgroun d 
turbulence itself can cause nonlinear damping to the waves ( Yan & Lazarianl 2002 ). 
Unlike hydrodynamical turbulence, MHD turbulence is anisotropic with eddies 
elongated along the magnetic field. T he anisotropy increases with the decrease of 
the scale ( IGoldreich & Sridharl 119951) . Because of the scale disparity, k\\ > k± ^> 
the nonlinear damping rate in MHD turbulence is less than the wave fre- 
quency k^v A , and it is g iven by ( Farmer & Goldreichl 2004 ; Yan & Lazarianl 2004 ; 
Beresnvak & Lazarian, 2008) 



(22) 

where L is the injection scale of background turbulence, and the k is set by the 
resonance condition k ~ k\\ ~ 1/ r£. 

There are various models for the diffusive shock acceleration. We consider here 
the escape-limited acceleration. In this model, particles are confined in the region 
near the shock where turbulence is generated. Once they propagate far upstream at 
a distance / from the shock front, where the self-generated turbulence by CRs fades 
away, the particles escape and the acceleration ceases. The characteristic length 
that particles penetrate into the upstream is D(p)/U. The maximum momentum 
is reached when D(p)/U ~ Z/41 Assuming / = KR s / t , where K < 1 is a numerical 
factor, one can get 

p max 3KA 2 UR sh 



mc \/l +A 2 vr 



(23) 



In particular, for A < 1 



mc 




i 2 



(24) 



7 The factor 1/4 arises from the following reason. As pointed out by Ostrowski & Schlickeiserl 
( 1996), the spectrum is steepened for small 1, i.e., IU /D(p) ~< 4 



Cosmic Ray transport in turbulent magnetic field 



21 



In the limit of low shock velocity, 



v A < U <C c 



KRsh 
c J 2LC cr at, 



1/4 



(25) 



we get 



Pmax 



3 ,2 



r c z v\ 



(26) 



At the Sedov phase (t > t sed = 250{E 5l / \n Ug)) l / 3 yr), where E51 = E SN /I0 5l eig 
and U9 = f/, /10 9 cm/s are the total energy of explosion and the initial shock velocity. 
In FigfSJ we plot the evolution of p max / (mc) during the Sedov phase. The solid line 
represents the results from Eqs.(l24V As we see, at earlier epoch when advection 
and streaming instability are both important, the evolution of p max does not follow 
a power law. For comparison, we also put a power law evolution in the same figure 
as depicted by E q.d26]i (dashed line). Our result is also larger than that obtained by 
IPtuskin & Zirakashvilil d2005l) since the wave dissipation rate is overestimated in 
their treatment. 




Fig. 5 Left : The energy of CRs that are released at th e shock at time t in the Sedov phase. Ou r result 
shows that the often assumed power law solution (see lGabici et al., 2009; Ohira et al.Ll20l"bT) is only 
realized in asymptotic regime as described in Eqs.l l25l2"6V It is also larger than the earlier result 
(dotted line) in Ptuskin & Zirakashvili (2005) where the damping of the waves by background 
turbulence is overestimated. Right: The spectrum of CRs at a distance r = 12pc after 1800 (solid 
line), 6000 (dotted line), 12000 (dashdo t line), 50000 ye ars (cross line). The Galactic mean is 
plotted as a reference (dashed line). From lYan et al.l d2012h . 



6.2 Enhanced scattering and streaming instability near SNRs 



The result from lYan et al.1 (12012b show that the local scattering of CRs has to be 
enhanced by an order of magnitude % = 0.05 in order to produce the amount of 
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7 ray emission observed. A natural way to increase the scattering rate is through 
the streaming instability. The enhanced flux of the CRs are demonstrated to gen- 
erate stron g enough instabil ity to overcome nonlinear damping by the background 
turbulence dYan et al.1 120121) . The growth rate in the linear regime is 



N(>E) 



VA 



(27) 



where v., is the streaming speed of CRs. The growth rate should overcome the damp- 
ing rate (eql22t for the instability to operate. The condition r gr > Q leads to 



v s > v A \ 1 + 



nv A 



N£2 , 



(28) 



The spatial diffusion coefficient adopted here, D rj v s L = xDism, satisfies this 
requirement. The growth and damping rates are compared in Fig\6\right. We see that 
the streaming instability works in the energy range needed to produce the observed 
y ray emission, proving that our results are self-consistent. 

Note that the case w e consider here is diff erent from the general interstel - 
lar medium discussed in Yan & Lazarianl d2004 and iFarmer & Goldreichl d2004 . 
namely, the local cosmic ray flux near SNRs is much enhanced (see Fig|5] right). 
Consequently, the growth rate of the streaming instability becomes high enough to 
overcome the damping rate by the preexisting turbulence in the considered energy 
range. 




Fig. 6 Left : The spectrum of Gamma ray emission from W28. The Fermi data are shown as dotted 
points jAbdo et alll2010t) . and the H.E.S.S. data are plotted as 'x' points jAharonian et al 
with error bars. Solid line is our result. Right. The growth and nonlinear damping rates of streaming 
instability. With the locally enhanced flux, the growth rate of streaming instability becomes much 
larger than the mean Galactic value so that it can overcome the nonlinear damping by turbulence 
for a wide energy range. This is consistent with our earlier treatment in which streami ng instability 
plays an essential role in the cosmic ray diffusion near SNRs. From lYanetaTlj2012h . 
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Table 3 Model parameters adopted 

a X r\ k \ a 

0.1~0.3 ~0.05 ~0.3 0.04 ~ 0.1 0.2 ~0.4 0.5 



7 Gyroresonance Instability of CRs in Compressible Turbulence 



(c) halo 



E(k) 




Fig. 7 Left: The spectral energy density of slab waves that is transferred from the large scale 
compressible turbulence via the gyroresonance instability of CRs. In the case that the instability 
grows up to the maximum energy rate allowed by the turbulence cascade, large scale turbulence 
is truncated at Xn and the wave amplitude E(k)dk ~ e'^ is given by Eq. OOt . Note that the picture 
is different from LB06, namely, the feedback on the large scale turbulence occurs only in some 
cases when the scattering is n ot sufficient to prevent the waves from growing to the maximum 
values jYan & Lazarian. 201 1]); right: CR scattering is dominated by compressible modes. For 
high energy CRs (>~10GeV), the scattering is due to direct interaction with fast modes; For low 
energy CRs, the interaction is mainly due to the gyroresonance instability induced by compression 
of magnetic fields. 



Until recently, test particle approximation was assumed in most of earlier stud- 
ies in which turbulence cascade is established from large scales and no feedback of 
CRs is included. This may not reflect the reality as we know the energ y of CRs is 
comparab le to that in turbulence and ma gnetic field (see iKulsrudi 120051) . It was sug- 
gested by lLazarian & Beresnvakl (120061) that the gyroresonance instability of CRs 
can drain energy from the large scale turbulence and cause instability on small 
scales by the turbulence compression induced anisotropy on CRs (see Fig(7]^e/0- 
A nd the wave generated o n the scales, in turn, provides additional scattering to CRs. 
In lYan & Lazarianl(1201 lb . quantitative studies was provided based on the nonlinear 
theory of the growth of the instability and the feedback of the growing waves on the 
distributions of CRs. 

In the presence of background compressible turbulence, the CR distribution is 
bound to be anisotropic because of the conservation of the first adiabatic invariant 
)j, = v\jB. Such anisotropic distribution is subjected to various instabilities. Waves 
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are generated through the instabilities, enhancing the scattering rates of the parti- 
cles, their distribution will be relaxed to the state of marginal state of instability 
even in the collisionless environment. While the hydrodynamic instability requires 
certain threshold, the kinetic instability can grow very fast with small deviations 
from isotropy. Here, we focus on the gyroresonance instability. Both the qualita- 
tive and quantitative studies in lYan & Lazarianl (1201 II) show that the isotropization 
rate is roughly ~ where r cr ,En are the instability growth rate and the 

wave energy normalized by magnetic energy, respectively. Pcr is the ratio of CR 
pressure to magnetic pressure, A is the degree of anisotropy of the CR momentum 
distribution. 

By balancing the rate of decrease in anisotropy due to scattering and the growth 
due to compression, one can get 



En 



r„ r vA 



Xcr = r p /e N . 



(29) 



where \>a is the Alfven speed, co,8v are the wave frequency and amplitude at the 
scale that ef fectively compre s ses th e magnetic field and create anisotropy in CRs' 
distribution dYan & Lazarianl 1201 lb . Since the growth rate decreases with energy, 
the instability only operates for low energy C Rs (~< lOOGeV , see F ig[71 right) due 
to the damping by the preexisting turbulence ( Yan & Lazarian . 201 ll) . 



7.1 bottle-neck for the growth of the instability and feedback on 
turbulence 

The creation of the slab waves through the CR resonant instability is another channel 
to drain the energy of large scale turbulence. This process, on one hand, can damp 
the turbulence. On the other hand, it means that the growth rate is limited by the 
turbulence cascade. The energy growth rate cannot be larger than the turbulence 
energy cascading rate, which is 1 /IpV^/vA/L for fast modes in low j3 medium and 
pv\/U for slow modes in high j3 medium. This places a constraint on the growth, 
thus the upper limit of wave energy is given by 



_[M\L i /{IA)f-\ J3<1 
£N -U/(l A A)f-\ P>1, (30) 

where 7 is the Lorentz factor and L; ~ 6.4 x 10~ 7 (fi/5^G)(10~ 10 cm 3 /n c ,)pc. The 
growth is induced by the compression at scales ^< Xcr- Therefore, in the case 
that r gr e reaches the energy cascading rate, fast modes are damped at the corre- 
sponding maximum turbulence pumping scale A«, = r p /EN (see Figf7]fe/f). If Xn, is 
larger than the original damping scale l c , then there is a feedback on the large scale 
compressible turbulence. This shows that test particle approach is not adequate and 
feedback should be included in future simulations. 
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8 Summary 

In this paper, we reviewed recent development on cosmic ray transport theories 
based on modern understanding of MHD turbulence. The main conclusion is: 

• Compressible fast modes are most important for CR scattering. CR transport 
therefore varies from place to place. 

• Nonlinear mirror interaction is essential for pitch angle scattering ( including 90 
degree). 

• Cross field transport is diffusion on large scales and super- diffusive on small 
scales. 

• Subdiffusion does not happen in 3D turbulence . 

• Our results are tested using input from turbulence simulations. 

• Self - generated waves are subject to damping by preexisting turbulence 

• Small scale waves are generated in compressible turbulence by gyroresonance 
instability, dominating the scattering of low energy CRs (< lOOGeV). Feedback 
of CRs on turbulence need to be included in future simulations. 
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